The effect of focusing on polarization qubits 
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When a photon with well-defined polarization and momentum passes through a focusing device, 
these properties are no longer well defined. Their loss is captured by describing polarization by a 
3x3 effective density matrix. Here we show that the effective density matrix corresponds to the 
actual photodetection model and we provide a simple formula to calculate it in terms of classical 
fields. Moreover, we explore several possible experimental consequences of the "longitudinal" term: 
limits on single-photon detection efficiency, polarization-dependent atomic transitions rates and the 
implications on quantum information processing. 
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I. INTRODUCTION 

Single-photon manipulations have recently become an 
integral part of quantum optics and they play an im- 
portant role in experimental quantum information pro- 
cessing 0. Already with current technology the single- 
photon states can be produced and manipulated quite 
reliably 0, and the efficiency of their production and 
detection is expected to rise in the future. Usually these 
states are considered to be eigenstates of momentum, and 
often the very notion of photon is synonymous with an 
elementary excitation of the electromagnetic field with 
a well-defined momentum and polarization. A single- 



photon state |k, a) 



vac) is an excitation of the 



plane wave mode of momentum k and polarization a and 
is distributed over the entire space. An approximately lo- 
calized photon PI is described instead by a superposition 



(1) 



operator- valued measures 0, 0] . Both experemental re- 
sults and theoretical calculations for classical fields and 
for quantum coherent states 0, 0, 0| establish new ef- 
fects that follow from the presence of a significant longi- 
tudinal electric field. This is often created by using fo- 
cusing devices, which are also an inevitable part of usual 
optical experiments. 

In this paper we apply the effective density matrix for- 
malism to quantum states that correspond to classical 
modes with a significant longitudinal component and dis- 
cuss its connection with experiment. The paper is orga- 
nized as follows. Sec. II reviews effective 3x3 density 
matrices and describes some of their applications. Sec. Ill 
presents a general formula that expresses an effective den- 
sity matrix in terms of classical modes. It is illustrated 
by an elementary discussion of the lens action on one- 
photon states. Sec. IV discusses a connection between 
the formal construction of effective density matrix and a 
simple photodetection model. Finally, Sec. V introduces 
possible experimental consequences of the "longitudinal" 
term in effective density matrices. 



where a denotes the helicity, the normalization is given 
by f dfi{k)\fcr{'k)\'^ = 1, while we adopt a non- 
relativistic measure (i/i(k) — (i^k/(27r)^. More generally, 
localized states may be described as mixtures of such 
terms. 

A typical spread in momentum is often very small, and 
polarization is approximately constant. As a result, po- 
larization is usually described by a 2 x 2 reduced density 
matrix which is formally equivalent to that of a qubit, the 
ideal unit of quantum information . To ensure validity 
of this approximation and to calculate possible correc- 
tions, a general notion of polarization density matrix is 
required. However, the standard definition of reduced 
density matrix fails for photon polarization 5], and it 
is possible to define only an effective 3x3 density ma- 
trix which corresponds to a restricted class of positive 
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II. EFFECTIVE DENSITY MATRIX AND ITS 
APPLICATIONS 

If one is interested only in polarization degrees of free- 
dom, it is tempting to define a reduced 2x2 density 
matrix by 



Pcr,a' 



rf/i(k)/,(fc)/;,(fc). 



(2) 



However, helicity eigenstates are defined only with re- 
spect to a given momentum. Intuitively, the polarization 
vectors for different momenta lie in different planes and 
cannot be superimposed. Under rotations each compo- 
nent acquires a momentum-dependent phase and hence 
the density matrix |(2Jl has no definite transformation 
properties This makes a standard density matrix a 

useless concept even when a fixed reference frame is con- 
sidered, since any POVM (positive operator-valued mea- 
sure 0,01) that describes an experimental setup must 



2 



have definite transformation properties at least under or- 
dinary rotations. Analysis of the one-photon scattering 
[T^ l gives another angle on the failure of this concept. 

The 3x3 matrix with the right transformation prop- 
erties may be introduced with the help of polarization 
3- vectors Q. A polarization state |Q!(k)) corresponds to 
the geometrical 3- vector 



a(k) = a+(k)e|^ + a^{k)e^ , 



(3) 



where ja+p + la- 1^ = 1, and the vectors correspond to 
the right and left circular polarization and satisfy k-ej = 
0. In this notation a generic one-photon state can be 
written as 



1*)= J d/.(k)/(k)|k,a(k)), 



(4) 



and the effective 3x3 description takes the form 

P^n = y"dMk)|/(k)|'a„(k)a„(k)*, (5) 

where m,n — x,y,z and the vector Q:(k) is given by 
Eq. (jSJ. Under rotations of the coordinate system this 
density matrix has a simple transformation law, 



p RpR^, 



(6) 



where R is the rotation matrix. The effective den- 
sity matrix p can be obtained with the standard state- 
reconstruction techniques from a family of POVMs 0. 
For example, its diagonal elements pa are the expectation 
values of the elements of the "momentum-independent" 
polarization POVM that consists of three positive oper- 
ators Ei that sum up to the identity, Ei = 1, 



(7) 



We discuss a relation between this POVM and the stan- 
dard photodetection model in Sec. IV. 

Let us now examine a wave packet which describes a 
nearly plane and nearly monochromatic wave. In this 
limit, /(k) of Eq. is strongly localized around some 
central value kg. By an appropriate transformation of 
the form ^ the effective 3x3 reduced density matrix 
can be put into a block diagonal form. Then the density 
matrix will have a 2 x 2 block with almost unit trace, 
which corresponds to the standard 2x2 reduced density 
matrix that would describe the state for /(k) cx 5(k— ko). 
As an example, consider the following wavepackets that 
are formed by the helicity eigenstates. 



l*±) = 



dMk)/(k)|k,e±), 



(8) 



where /(k) satisfies the above criteria. Calculating the 
effective reduced density matrix we have 
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(9) 



where f2 4; 1 is roughly the ratio of a typical width of 
the wave packet to jkol, and p_ = p+. Its exact form de- 
pends on the detailed shape of /(k). Thus we arrive to 
the conclusion that integrating out the photon's momen- 
tum leads to polarization states that are neither pure nor 
perfectly distinguishable j^. There is a non-zero proba- 
bility to identify a state as a p_ state and vice versa. 
The probability of making such an error given a perfect 
equipment and measurement scheme [T3 | will be denoted 
by Pe and is given by 



PEip+,P-) 



itr|p+ - p_ 



2" 



(10) 



Although conceptually the implications of the above con- 
clusions are profound, when considering a wave packet 
with a very narrow distribution in momentum /(k), 
the parameter f2 is very small. Consider an electro- 
magnetic beam propagating along the 2-axis and with 
the Gaussian distribution in intensity in the (xy) plane 
that is given by /(r) = /q exp(— r^/r^). Since the dis- 
tribution in momentum /(k) is essentially the Fourier 
transform of the electric field, we expect that the radial 
spread in momentum will also be of the Gaussian form 
/(k) oc /i(fc^)exp(-fc2/2A2) where - 1/t, and /i is 
some function of k^- Assuming also the Gaussian distri- 
bution in wavelength fi{kz) ~ exp(— (fcz — fco)^/2A^) we 
have 

/(k) = 7Vexp(-(fc, - A;o)V2A2)exp(-A;2/2A2). (11) 

For this distribution in momentum, the parameter f2 is 
given by = Ar/ko + 0{A'^/k^). Taking r to be of 
the orders of 10~'^m, and a wavelength of 5 x 10~''m, we 
get ri ~ 5 X 10~^, thus rendering the effect negligible. 

However, the momentum spread becomes substantial 
when a beam undergoes focusing. A classical electromag- 
netic plane wave that passes through a converging lens is 
no longer plane or transversal 0, • A substantial lon- 
gitudinal field component is present as well as the effect 
described above. 



III. THE EFFECTS OF FOCUSING 

In the analysis of the influence of a lens on quantum 
states it is important to bear in mind that it is incapable 
of turning a pure state into a mixture unless informa- 
tion is lost. Indeed, pure incoming state l^'in) trans- 
forms into a pure outcoming state l^'out)- The same 
information is encoded differently in these two states, 
and the effective density matrix captures our inability 
to access all of it. The lens action is analyzed by mode 
matching. To achieve this, we find a mode decomposi- 
tion of solutions of the corresponding classical equation 
(in the case of photons these are the electromagnetic 
wave equations for the vector potential. Throughout 
this paper we use the Coulomb gauge). With each mode 



(27r)«6^e- 



-i(ujt k x) ^ creation operator al^_^_ is as- 



sociated. The resulting quantum state |k, ±) 



vac) 
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is normalized as (k'cr'|k, cr) = S'^^^'k — k')Sca' ■ AH clas- 
sical solutions are of the form 



A(t,x) = A{t,x)a{t,x) = J dn{k)a{k)A{k)e~ 



■2(cj£— k-x) 



(12) 



with E = —A. The frequency uj = \k\c, each field compo- 
nent is split into the field strength A{k) and the transver- 
sal polarization part Q:(k), k • cy.{k) = 0, |a(k)| = 1. The 
corresponding norm is 



ll^ll 



d^lik)\Aik)\ 



(13) 



Therefore, a normalized positive energy solution of the 
form 1)12(1 corresponds to a one-particle state 

I*) = E / rf/^(k)a.(k)/(k)4jvac) 

dM(k)/(k)|k,a(k)}, (14) 



where the coefficients aa are defined by Eq. 0, and 
f(k) = ^(k)/||yl||. The transformation |*out) is 

obtained from the transformation of the incoming modes 
Aj^ into the outcoming modes A™* 0, Q • When the 
classical solution corresponding to the quantum state is 
known, the effective polarization density matrix Eq. (jSJ 
can be calculated as follows. We note that 



^„(x,t)A;,(x,i)(i3x= Ar,ik)AU^)dfiik), (15) 



for all t and n,m - 
can be written as 

P = 



The density matrix (jSJ, thus, 



/d^(k)A(k)At(k) 
/d/.(k)|^(k)P 



(16) 



We stress that this is a general expression independent 
of approximations that are generally used to calculate 
classical solutions. 

The analysis of the lens action should be performed in 
the vector diffraction theory and reliable estimates 
of a field in the focal regions follow the techniques of 
Richards and Wolf In some cases it is even possible 
to get exact solutions of Maxwell equations Then, us- 
ing Eq. H16() we find the outcoming quantum state whose 
effective reduced density matrix is calculated according 
to Eq. ©. 

To illustrate the importance of the longitudinal term, it 
is sufficient to consider an incoming state with a definite 
momentum and polarization |k, e^), i.e. to approximate 
the corresponding classical field by a plane monochro- 
matic wave. As we showed above, this is a good approx- 
imation for incoming wavepackets that are the actually 
used in experiments. Moreover, it is sufficient to use ray 
tracing, while more refined calculations should be used in 
conjunction with concrete experimental schemes. Hence 



we take the outcoming classical field as a spherical wave 
that converges to the geometric focus of a thin lens with 
the focal length /, as illustrated on Fig. 1. The polar- 
ization direction and the field strength at each point are 
calculated using the eikonal equation and ray tracing |l5j| . 

The set-up is schematically presented on Fig. 1, while 
calculations are given in Appendix A. For the incoming 
states |k, e^), in the leading order the outgoing states 
are 



P± 



74) I ±i 
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(17) 

This is the only approximation which is consistent with 




optical axis 



FIG. 1: Transformation of an incoming plane wave into a 
spherical wave by the focusing system. 

use of Gaussian optics for mode calculation. All three 
eigenvectors are non-zero (even if, of course, the one that 
is associated with an ideal circular polarization domi- 
nates). The probability of error in distinguishing between 
and p_ is proportional to the square of the numerical 
aperture. 



(18) 



Since 0m « /// where / is the aperture radius, for 
the moderate numerical aperture of the order of 10~^, 
the probability of error gets to the level of percents, 
compared with a case with no lens present. The ef- 
fects of "longitudinal" polarization will be even more pro- 
nounced when quantum states that correspond to classi- 
cal doughnut-shaped @ and other exotic modes are pro- 
duced. This will be enhanced by going to higher numer- 
ical aperture, where the longitudinal field can contain 
nearly 50% of the total beam power in the limit NA 1. 



IV. MEANING OF p 

While the POVMs that were described in 6] are le- 
gitimate theoretical constructions, the resulting 3x3 ta- 
ble becomes experimentally relevant if the POVM corre- 
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spends to some detection model. We establish this corre- 
spondence. Let us compare predictions of the POVM E^, 
Ey, Ez,J2 = 1 that gives the diagonal elements of the 
effective density matrix pjj = tr(£'j|<I>)(<I>|) '3, with 
the detection probabilities that can be obtained from 
the first-order perturbational calculations of the follow- 
ing model detector. The operators Ej are given explicitly 
in Appendix B. 

To facilitate the comparison we use an alternative form 
of the diagonal elements of the density matrix as in 
Eq. 



where 



■r 



dxdydz, W = '^Wj 



(19) 



(20) 



We discuss here semiclassical detection theory, since in 
the leading order the results of the full theory agree with 
the semiclassical one Q . In the latter the electromagnetic 
field is considered classically, while the photoelectrons are 
treated quantum-mechanically. The interaction term is 
given by p • A, where A is a classical vector potential 
(whose direction in the polarization gauge is given by the 
polarization vector a) and p is the electron's momentum 
operator. Accordingly, in the following we use a classical 
language to describe the electromagnetic field. The area 
5 of a planar detector is assumed to be much larger than 
the cross section of the beam. We model the detector's 
sensitivity to the wave's polarization by restricting the 
electron momentum to lie only along a chosen direction. 

Let us first consider a circularly polarized monochro- 
matic beam in the paraxial approximation |l7j |. Assum- 
ing that it propagates along the z-axis, we have 



E(x,t) « [E{x,y)e 



dE ^.dE 
dx dx 



(21) 



B « TiE, (22) 

where the basis polarization vectors are = 
^(1, ±z, 0) and the beam radius t is much larger than 
a typical wavelength, rfco S> 1. 

Let the planar detector absorb the field along the j-axis 
(j = x,y,z) and locate it at z = zq. Hence, electrons' 
excitation rate is proportional to J|A(x,?/,zo) • 
Assuming a finite detection time A (or a beam of finite 
duration A, with the fields that are given by the cor- 
responding superpositions of E and B having different 
frequencies with some weight function f = /(lu — loq)), 
the excitation probability is proportional to 



\A{x,y,zo) ■il'^dSdt. 



(23) 



Since the detector is planar, we can change the integra- 
tion variable from t to z we get 

/j = A y I A(a;, y, zo) ■ j\^dxdy (24) 

A{x,y, z) ■ dxdydz. (25) 



z=za 

I 

x,y J z—zq~ Ac 



Hence Ij — Wj , and the diagonal elements of p are indeed 
related to the photocurrent as 



Ps 



(26) 



Now consider a spherical wave which represents the 
EM field after the lens, as in Eqs. We 
again consider a planar wave detector and detection time 
A. This time, however, the normal to the detector plane 
is not parallel to the Pointing vector, and thus Ij/I 
are different from Wj/W. Nevertheless, as will be seen 
shortly, these ratios agree up to the order 9'^, while for 
the realistic values of the interesting effects (such as 
the predicted error probability, Eq. (|18(l ') are of the order 

Let us assume that the detector is located aX z = zq 
behind the focus. We use spherical coordinates with the 
origin in the focus and = at z = zq, and polar coordi- 
nates in the detector plane. In this plane cf) is the same as 
in the spherical coordinates and r = zotan^. From the 
result of Appendix A, the field strength E can be written 
as 



E = E{R,t 



1 



1 



(cos 0)3/2 r 



e{9) 
r 



- = (27) 



and since we are working with a well defined frequency, a 
similar decomposition is possible for A, i.e., A = a{9)/r. 
The point {9, cj)) on the detector plane is at a distance 
Zo/ cos 6* from the focus. The area element is 



dS = 



smf 



cossr" 

so the detection probability is proportional to 



(28) 



Ij^^J HO)\^\a{0,(l}) • j)ptan6'd6'd(/). (29) 
On the other hand, the integration over the shell Ac gives 



\a{9)\^\cx{9,cl>)-j)\^ sin 



The normalized detection probabilities 



(30) 



(31) 



differ therefore from the matrix elements of pjj = Wj/W 
by the presence of additional factors 1 / cos 9 in each of 
the integrals. However, when the results are expanded in 
terms of 9^, the difference between these two expressions 
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is only of the order of or higher. Accordingly, the lead- 
ing order expansion for the probability of error, Eq. (|18|) . 
remains the same for the model detection scheme that 
was described above. This discrepancy highlights the fact 
that different detection procedures lead to different polar- 
ization density matrices |12| |. Moreover, with a slight 
abuse of the language, we note that the resulting effec- 
tive mixing is caused by the shape mismatch between the 
wave front and the detector. 



V. SUMMARY AND OUTLOOK 

We have shown that an effective 3x3 polarization 
density matrix, previously introduced on formal grounds, 
indeed has a direct experimental significance. It should 
be noted, however, since there is no general polarization 
density matrix, different detection procedures may lead 
to different effective constructions. 

Presence of a significant longitudinal part in the effec- 
tive density matrix imposes limits on detection efficiency. 
The interaction between single atoms and electromag- 
netic field, either classical or in a quantum coherent state, 
is affected by focusing ^ . The structure of the ef- 
fective density matrix shows that this will be the case 
also for single photon states. Polarization-sensitive tran- 
sitions in the atoms will be suppressed by strong focusing. 
For low numerical apertures we expect these effects be 
proportional to its square, reaching saturation on higher 
levels when N A — > 1. From the point of view of quantum 
information theory, the fact that these density matrices 
are inevitably mixed, actually implies that polarization 
qubits are always noisy. This intrinsic noise should be 
taken into account in the analysis of the physical realiza- 
tions of quantum computing and in the security analysis 
of quantum cryptographic protocols. Similarly, sensitiv- 
ity of photon-atom interaction to the focusing will affect 
the efficiency of trapped-atom based quantum memory 

Q 
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APPENDIX A 

Assume that the optical axis is the z-axis, and approx- 
imate the incoming classical wave as a plane wave with 



k = z and polarization a. Two systems of coordinates 
will be used in the following. The plane z = is the 
plane where the lens is situated. In it we define polar 
coordinates (r, cj)) , where r is measured from the optical 
axis in the 2 = plane. The unit vectors f , <p will always 
lie in the z — plane. The spherical coordinates (i?, 9, cj)) 
are calculated from the focus (the angles cj) in both coor- 
dinate systems are the same). The intersection between 
the optical axis and z = plane has r = f and 6 = 
coordinates, while the lens aperture is bounded by 0m • 
The relation 

r^/tane, (32) 

holds on the plane z = 0, while ^m, 

tan^m^Z//, (33) 

is determined by the nominal focal ratio jlSI] F = //2Z, 
where I is the radius of the entrance pupil. In the Gaus- 
sian approximation 0m equals to the numerical aperture. 
Assuming that after passing through the lens the field 
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FIG. 2: Changes in a linear polarization according to the 
ray-tracing model. The direction and polarization labels k 
and a of the ray are transformed into K{6,(f>) and a.'(9,cf>), 
respectively. 



is that of a perfect spherical wave, the ray that hits the 
z = plane at (r, cj)) is deflected to the direction 

k'(6') = — sin 6*? -l-cos 6'z = — sin0(cos 0x-|-sin 0y)-f cos 6'z, 

with 9 given by Eq. H32|l . In passing through a homoge- 
nious dielectric medium polarization directions are paral- 
lel transported along each ray. At the boundaries the di- 
rection that is parallel to the incidence plane is refracted 
to remain transversal to the wave vector, while the per- 
pendicular component is unchanged [isl ll^ . The inci- 
dence plane is spanned by the vectors z and r, so that 
the parallel and perpendicular components of polariza- 
tion are f and 0, respectively. We decompose the polar- 
ization of the ray that passes through the z = plane at 
(r, 4>), in terms of f and <j). For linear x-polarization it is 
given by 

eta; = X = cos 0f — sin (l>4>. (35) 
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Having passed through the z = plane, the new direction 
of polarization for the above ray becomes 



— cos 0(cos 6r + sin 9z) — sin c 



(36) 



were again 9 is related to r by Eq. H32|l . The same cal- 
culation can be carried out for the linear y polarization, 
giving 



a. 



— sin 0(cos Or + sin 9z) + cos ( 



(37) 



It is easy to see that right and left circular polarizations 
are preserved up to a phase, e^^'-^^'^-'e^, , where the 

precise form of g is irrelevant. Fig. 2 illustrates these 
changes. 

To complete the classical description of the field as in 
Eq. (|12|l . we need the field strength i?(k). Calculations of 
the intensity are based on the intensity law of geometrical 
optics 01, 



(38) 



where E, dS and E', dS' are the field strength and the 
area element at the respective wavefronts. Taking the 
initial field strength to be unity, and considering that 
the wavefronts before the lens are planar. 



dS = 27rr(e)dr(e) = 2T:^^fde, 

and that after the lens they are spherical, 
dS' = 2T:smeR^de, 

we get 



E\R,e) = 



1 



/ 

(cosf?)3/2 r 



(39) 



(40) 



(41) 



Since both the plane wave and the spherical wave are 
non-normalizable, we obtain the density matrices from 
the two-dimensional integration over the angular parts 
of the volume integral. For the incoming states |k, e^), 
taking into account that w = |k'|c = const, the outgoing 
states are 



_ C=oIe% sin 0ded^e±i9, ^)eii9, ^)/ cos^ 9 
^~ J dedijjsind/cos^e 

whose explicit form is given by Eq. H17|l . 



(42) 



of our construction. The POVM itself is build only with 
the physical polarization states. Allowing for longitudi- 
nal polarization makes it possible to define a polarization 
state along an arbitrary direction, say the x-axis, as 

\k) ^ x+{k)\e+) + x^{k)\e^) + x,{k)\ei), (43) 
where x±{'k) = • x, and x^(k) = x • k. Note that 
(x|y) = X ■ y = 0, whence 



|x)(x| + |y)(y| + |z)(z| = l. 



(44) 



A projection operator that corresponds to the direction 
X is 



= ® Ip = ® J dAi(k)|k)(k|, 



(45) 



where Ip is the unit operator in momentum space. The 
action of on a physical state j'l') follows from Eq. H43|) 
and (e^ kk) ~ ^- Only the transversal part of |x) appears 
in the expectation value: 



(vI/|P,|vI/) = J dMk)|/(k)n:r+(k)a;(k)+x_(k)a*_(k)p. 

(46) 

Define the transversal part of |x): 

|k,b.(k))EE(|e+)(e+| + |e^)(e^|)|x) 

-x+(k)|e+)+x_(k)|e^), (47) 

and likewise |bj,(k)) and |b2(k)). These three vectors are 
neither of unit length nor mutually orthogonal. 

Finally, a POVM element E^ which is the physical 
part of Px , namely is equivalent to P^ for physical states 
(without longitudinal photons) is 



Ex= J dM(k)|k,b,(k))(k,b,(k)|, 



(48) 



and likewise for other directions. The operators E^, Ey 
and Ez indeed form a POVM in the space of physical 
states, owing to Eq. ll^ . It then follows from Eq. jT7|) 
and similar definitions for the other directions that, for 
any k. 



APPENDIX B |b.(k))(b.(k)| + |b,(k))(b,(k)| + |b,(k))(b,(k)| = l^k, 

(49) 

The required POVM is introduced as follows. The Ion- where Ij^k is the identity operator in the subspace of 
gitudinal photons are used to define the necessary steps polarizations orthogonal to k. 
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